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ABSTRACT 

Recent observations of the compact source embedded within the supernova remnant RCW 
103 rekindle interest in the origin of this object's emission. We contrast several models in which 
neutron-star cooling powers RCW 103. Specifically, either the presence of an accreted envelope 
or a sufficiently intense magnetic field can account for the X-ray emission from this object. 

Subject headings: stars: neutron — stars: magnetic fields — radiative transfer — X-rays: stars 



Introduction 



Soon after the X-ray source IE 161348-5055 was first detected by the Einstein observatory ( Tuohy fc 



Garmire 1980) near the center of the supernova remnant (SNR) RCW 103, Tuohy et al. (1983) proposed 
that this source is an isolated neutron star emitting thermal radiation. Optical and radio observations 
have failed to identify a counterpart ( Tuohy fc Garmire 198C| ; Tuohy et al. 1982; Dickel et al. 1996; Kaspi 



et al. 1996), bolstering the interpretation of this source as an isolated neutron star. Subsequent X-ray 
observations with Einstein and i? 05*^4 T" have not all confirmed the initial detection (Tuohy & Garmire 198C| ; 
Becker fc Triimpcr 199^ ). 



Using recent observations of IE 161348-5055 with the ASCA observatory and archival data from 
ROSAT, Gotthelf, Petre fc Hwang (1997)| verify the existence of this source and refocus attention on the 
interpretation of its emission. After subtracting a model for the emission of the surrounding SNR, [Gotthelf 



Petre fc Hwang (1997)| find that the point source spectrum is well described by a blackbody having a 



characteristic temperature kT — 0.6 keV and a flux of 6 x 10~ erg s~ cm~ . Estimates of the distance to 
RCW 103 vary from 3.3 kpc ( [Caswell et al. 1975 ) to 6.6 kpc (Leibowitz fc Danziger 1983). Combining these 
values yields an estimated luminosity of 8d^ 3 x 10'^^ erg s~^ and an effective emitting area of 7c?| 3 x 10^'^ 
cm^ where c?3.3 is the ratio of the true distance to the X-ray source to 3.3 kpc. This is less than a percent of 
the total surface area of a neutron star. So, unless the emission originates from a tiny hotspot, a blackbody 



cannot account for the emergent spectrum. Gotthelf, Petre fc Hwang (1997) also find no periodic variation 



in the flux greater than 13 % of the mean count rate. Variation of this order or larger would be expected 
from a rotating neutron star emitting from a small portion of its surface unless the hot spot coincides with 
the rotation axis, the object's period is outside the range explored, or gravitational defocusing smooths the 
periodic signal. 



Several models for this object have been proposed since its discovery. [Gotthelf, Petre fc Hwang (1997) 



argue that the object is unlikely to be a cooling neutron star, a plerion, or a neutron star with an ordinary 



^Current address: Theoretical Astrophysics, mail code 130-33, California Institute of Technology, Pasadena, CA 91125 



- 2 - 



companion. The dismissal of these models prompted Popov (1997) to argue that the emission from 
IE 161348-5055 is powered by accretion onto a neutron star in a binary with another compact object. 
Unless the magnetic field of the neutron star is exceptionally weak {B < 10^ G), it will channel the accreted 
material onto the polar caps producing hotspots and variability. For an apparently young object to have 
such a weak field, he argues that the neutron star is not coeval with the remnant, but that the remnant 
resulted from the supernova of the binary companion. 

In this Letter, we revisit models of IE 161348-5055 which account for its emission through neutron star 
cooling. In the first, a neutron star cooling through an accreted envelope naturally results in a spectrum 
which greatly departs from a blackbody. The second model, an ultramagnetized cooling neutron star, 
results in anisotropic emission from a hotspot with a spectrum which qualitatively resembles a blackbody. 



2. Analytic Models 

The insulating envelope of a neutron star may be modeled analytically if the magnetic field is sufficiently 



strong or weak. We use the models of Hernquist & Applegate (1984) to describe the heat transport through 



an unmagnetized hydrogen envelope. The ultramagnetized models of Heyl & Hernquist (1997a) describe the 



relationship between core temperature and transmitted flux for ultramagnetized envelopes {B > 10^^ G). 



For such intense magnetic fields, these analytic models agree well with fully numerical calculations (Heyl 



fc Hernquist 1998|). Both sets of calculations adopt a plane parallel approximation to solve the thermal 



structure equation in the envelope and assume that the passage from the non-degenerate to the degenerate 
regime is abrupt. 

The equations are solved from the surface toward the core using a zero flux boundary condition 



( Schwarzschild 1965 ). The core temperature is a function of the combination F/gs and ip (in the magnetized 
case). Here, F is the transmitted heat flux, is the surface gravity, and -0 is the angle between the radial 
and field directions, and all of these values are taken to be in the frame of the neutron star surface. In the 
zero field limit, the models are independent of angle and in the ultramagnetized limit F/gs oc cos^ ip for a 
fixed core temperature. 



2.1. An Accreted Envelope 



As a baseline model, we consider a neutron star that has accreted sufficient unprocessed material 
since its birth to have a light-element envelope. We consider a hydrogen layer which extends all the way 
down to a density p ~ 10^" g cm~'^. The total mass of accreted hydrogen would be ~ 10~*Mq ( Heyl & 



Hernquist 1997bD 



To obtain the core temperature as a function of transmitted flux, we recalculate the models of 



Hernquist fc Applegate (1984) with the atomic number {Z) and weight (A) equal to unity. Since Hernquist 
fc Applegate (1984) examine iron envelopes for which Z + 1 « Z, we have to make several alterations to 
their results to account for the pressure contributed by the ions in the non- degenerate regime. 



Specifically, if the conductivity has a power-law form. 



Kq 



(1) 
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then the relation of temperature to density in the non-degenerate regime follows a solution such that the 
conductivity is a constant 

_a + f3F{Z + l)k 



a gs Anip 
The relationship between T and p takes the form, 

+ 13 F {Z + l)k 1 



T = 



a gs 



Kq 



V/3 



(2) 



(3) 



We assume that free-free scattering dominates the opacity through the non-degenerate portion of the 
envelope; consequently, 

„ ^ 13 16(7 196.5^2 g 

Since electrons dominate the pressure in the degenerate regime, we use equations (2.22) through 
(2.29) of Hcrnquist & Applegate (1984) without alteration to determine the core temperature given the 
temperature at the onset of degeneracy calculated using Equation ^. The more detailed unmagnetized 
models of Potekhin, Chabrier fc Yakovlev (1997) show that neutron stars with partially accreted envelopes 
exhibit cooling evolution between that of objects with fully accreted envelopes and standard cooling 
scenarios. 



2.2. An Ultramagnetized Iron Envelope 



We use the calculations for ultramagnetized iron envelopes described in Heyl fc Hcrnquist (1997a) . As 
in the unmagnetized case, we assume that free-free scattering dominates the opacity in the non-degenerate 
regime. However, we apply an anisotropy factor to account for the effect of the magnetic field on the 
scattering rates. We use the results of Pavlov fc Panov (1976) and Silant'ev fc Yakovlev (1980) to estimate 
this effect. 

In the degenerate regime, we proceed analytically to a density where the first Landau level fills. Below 
this density, the relationship between chemical potential and density is analytically invertible in the fully 
degenerate limit; consequently, for an iron envelope with B = 10^^ G, we can analytically integrate the 
structure equations up to a density of 1.5 x 10^° g cm"'^ using the conductivities for the liquid and solid 



phases presented by Hcrnquist (1984) , For the magnetized envelopes the emission is anisotropic along the 
surface. The average flux over the entire surface is a factor of 0.4765 times its peak value at the magnetic 



poles (Heyl fc Hcrnquist 1997a) 



With each of these models we determine the relationship between Toff and Tc for lO^K < Tcff < 10''K 
and calculate the cooling curves as described in Heyl fc Hcrnquist (1997c) . Figure |^ traces the cooling 
evolution for both ultramagnetized and unmagnetized iron envelopes, and unmagnetized hydrogen envelopes. 
We find that an intense magnetic field increases the emitted flux from an neutron star at a given time 
during the neutrino-cooling epoch by an amount sufficient to account for the observed luminosity of the 
point source in RCW 103 of ^ 10'^'* erg s~^, if the neutron star is approximately 1,000 years old. This age 



estimate is consistent with the observations of the remnant ( Tuohy et al. 197E ; Nugent et al. 1984 ; [Carter 



Dickel fc Bomans 1997 ). However, the effective temperature, kT = 0.2 keV, falls short of the observed 



characteristic blackbody temperature of 0.6 keV. 
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The presence of a fully accreted envelope dramatically increases the effective temperature for a given 
core temperature. For a 1,000-year-old neutron star the effective tcmpcratiue kT — 0.3 kcV is only a factor 
of two below that observed. The total luminosity is ^ 10'^''' erg s^^, greatly exceeding that observed. 

Due to gravitational fractionation, a neutron star which has accreted any hydrogen will have a hydrogen 
atmosphere. The presence of a hydrogen atmosphere shifts the emission blueward from that of a blackbody. 



2.3. Analytic power-law atmosphere 

In the LTE structure and NLTE spectrum formation limit of studying an atmosphere, it is 
straightforward to derive the spectrum for power-law conductivities. In general the conductivity and its 
associated opacities are given by 

^ = «o-, K=—-^ndKE = fhh ) (5) 

where Re is the opacity as function of photon energy and /(7) is obtained by calculating the Rosseland 
mean 

dB / EV ,^ / dB ,^ 15 f°° x^e'^x^ 



where B{E) is the Planck function for the intensity of blackbody radiation. 

We know that along a solution to LTE structure equations dT/ dz = F/k where both F the flux and k 
are constant. The optical depth at a given energy is given by 

riE.T, = I' -M. = £'.4^r = £ f = /(.) jl-H (Z.)' (I^f 

where we have changed variables, assuming the radiative zero solution, and substituted for ke and 
F = ctT^j.. Wc define Te = E/k. 

In the limit of a blackbody atmosphere, we have 7 = 0, 7(7) = 1, and T = Teg where the spectrum 
forms, therefore we take r = 4/3 to determine the temperature at which photons of a given energy have 
their effective photosphere. After some rearrangement we obtain, 

T.t. . 4 X-V(4+.) x./(4+.) 

7-=i/W4T^j [tTJ ■ 

To calculate the spectrum, we make the following three assumptions 

1. The Rosseland mean opacity is unaffected even though the spectrum diverges from a blackbody 
through the atmosphere, 

2. If the optical depth for a photon of a given energy to escape to infinity is less than 4/3, it is completely 
free, and 

3. If the optical depth is greater than or equal to 4/3, the photon is drawn from a blackbody distribution 
at the appropriate temperature. 
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At each layer of the atmosphere, flux conservation is imposed by scahng the blackbody contribution to the 
flux such that the total flux is equal to crT^g- With this prescription, it is straightforward to derive the 
emergent spectra as depicted in left panel of Figure |^. Even in this simple model, we find that an energy 
dependent opacity can shift the emergent spectra blueward of a blackbody. We find for 7 = 3, appropriate 
for free- free scattering ( Kippcnhahn fc Weigert 1990D that the peak is shifted blueward by nearly 40 %, 



which would result in an underestimate of the emission region by a factor of about four. 

To understand the dependence of the models on the first assumption, we can calculate the function 
7(7) using the emergent spectrum rather than a blackbody. By iterating from five to ten times, we find the 
appropriate values of f{j) to within one part per thousand. These iterated spectra are depicted in the right 
panel of Figure |[ For all values of 7 the blueward shift is less pronounced than in the direct prescription. 
For example for 7 = 3, the peak is shifted blueward by 25 %, which would result in an underestimate of the 
emitting area by a factor of 2.5. A self-consistent power-law atmosphere would require a recalculation of 
the mean opacity at each depth of the atmosphere. Within the model, according to assumptions (2) and 
(3), the spectrum at a given depth consists partly of a blackbody and partly of the emergent spectrum; 
therefore, we expect that a spectrum calculated self-consistently throughout the atmosphere would fall 
between these two limiting cases. 

More detailed modeling of hydrogen atmospheres supports our conclusions. As shown in Figure ^, 



this effect is even more pronounced when one examines the detailed calculations of Zavlin, Pavlov 



fe Shibanov (1996). They find that if one fits a blackbody to the emergent hydrogen spectrum one 
will overestimate the effective temperature by an even larger factor, ~ 4.2, than found in our models. 
Consequently, an ultramagnetized neutron star with an iron envelope and a hydrogen atmosphere can 
account for both the observed flux and characteristic temperature of IE 161348-5055. 



Conclusions 



We find that a young neutron star cooling through a strongly magnetized or a partially accreted 



envelope can account for the observed emission from IE 161348-5055. The detailed models of Potekhin 



Chabrier fc Yakovlev (1997) support the conclusions that we have found analytically in this Letter. 
Furthermore, the estimates of the emitting area of IE 161348-5055 support the conclusion that the spectrum 
from this object is significantly harder than a blackbody and possibly results from emission through a 
hydrogen atmosphere. 



The recent discovery by Torii et al. (1998) of a 69-ms X-ray pulsar (J161730-505505) in the vicinity of 



SNR RCW 103 complicates the evaluation of the possible models. It has a spin-down age of 8.1 x 10 yr. 



several times larger than that of the remnant. Torii et al. (1998) examine the possibility that the X-ray 
pulsar is associated with the remnant, and find that a kick velocity of 1300^3.3^8 i[ required 
to explain its current position relative to the center of the supernova remnant. Such a large supernova 
kick velocity is uncommon but has been observed for other pulsars (Lyne & Lorimer 1994). However, 



when this object is compared with other similar rotation-powered, plerionic pulsars, it is a factor of ten 



underluminous. We would argue with Gotthelf, Pctrc fc Hwang (1997) that this source is a heavily absorbed 
background object, possibly a rotation-powered, plerionic pulsar as Torii et al. (1998) suggest but located at 
a distance ~ 10 kpc. A reanalysis of the spectrum studied by Torii et al. (1998) may be able to determine 
the interstellar cohmm density to J161730-505505 and verify its status as a background source. 

I 1 

Gotthelf, Petre & Hwang (1997) failed to find flux variations at a level of 13 % over a wide range of 
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periods. Although the total flux from a magnetized neutron star may not vary at this level because of 



gravitational defocussing {e.g. Heyl & Hernquist 1997a), the magnetic field causes the atmospheric emission 



to be highly anisotropic (ghibanov et al. 1995; Rajagopal, Romani fc Miller 1997), so its apparent lack of 



variability may indicate that it is only weakly magnetized {B ^ 10^^ G) or simply that the geometry is not 
conducive to large flux variations. Observations of this object with AXAF should be able to distinguish 
between these models by determining the spectral shape of IE 161348-5055. 

We argue that the X-ray source in the supernova remnant RCW 103 is simply the natural end product 
of stellar evolution through a supernova: an isolated, cooling neutron star. 

We would like to thank the referee G. Chabrier for useful suggestions that improved the original 
manuscript. The work was supported in part by a National Science Foundation Graduate Research 
Fellowship and Cal Space grant CS-12-97. 
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Fig. 1. — The cooling evolution calculated analytically for iron and hydrogen envelopes. The results for 



unmagnetized iron envelopes are taken from Hernquist fc Applegatc (1984) . The upper cross-hatched region 
shows the fitted blackbody temperature of IE 161348-5055 and the acceptable range in age from 1,000 to 
3,000 years. The lower shaded region depicts the luminosity of the object estimated from observations. 
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Fig. 2. — Emergent spectra from the power-law atmospheres. The left panel shows the direct calculation, 
and the right depicts the iterated results. The solid line traces a blackbody spectrum, the dotted line is for 
a 7 = 1 opacity, the short-dashed line is 7 = 2, and the long-dashed line is 7 = 3. The heavy solid line shows 
results for a hydrogen atmosphere with T^s = 10^'^ K and gs = 2.43 x 10^'* cm s^^ from Zavlin, Pavlov fc 



Shibanov (1996) 



